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ABSTRACT 



We present the X-ray source detection procedure that we 
have developed for the purpose of assembling and charac- 
terizing controlled samples of cluster of galaxies for the 
XMM Large Scale Structure Survey. We describe how 
we model the selection function by means of simulations: 
this leads us to define source classes rather than flux lim- 
ited samples. Focussing on the CFHTLS Dl area, our 
compilation suggests a cluster density higher than previ- 
ously determined from the deep ROSAT surveys above 
a flux of 2 x 10 -14 ergcm -2 s~ 1 , We also present the 
L-T relation for the 9 brightest objects in the area. The 
slope is in good agreement with the local correlation. The 
relation shows luminosity enhancement for some of the 
0.15 < z < 0.35 objects having 1 < T < 2 keV, a 
population that the XMM-LSS is for the first time sys- 
tematically unveiling. 
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1. INTRODUCTION 



The question of cosmic structure formation is substan- 
tially more complicated than the study of the spherical 
collapse of a pure dark matter perturbation in an expand- 
ing Universe. While it is theoretically possible to pre- 
dict how the shape of the inflationary fluctuation spec- 
trum evolves until recombination, we hardly understand 
the subsequent galaxy, AGN and cluster formation be- 
cause of the problems of non-linear growth and the feed- 
back from star formation. Consequently we cannot use 
the statistics of visible matter fluctuations to constrain the 
nature of the Dark Matter and Dark Energy without devel- 
oping this understanding of non-gravitational processes. 
Clusters, the most massive entities of the Universe, are 



a crucial link in the chain of understanding. They lie 
at the nodes of the cosmic network, and have virialized 
cores, but are still growing by accretion along filaments. 
The rate at which clusters form, and the evolution of their 
space distribution, depends strongly on the shape and nor- 
malization of the initial power spect rum, as well as on th e 
Dark Energy equation of state (e.g. Rapetti et al. (2005)). 
Consequently, both a 3D mapping of the cluster distribu- 
tion and an evolutionary model relating cluster observ- 
ables to cluster masses and shapes (predicted by theory 
for the average cluster population) are needed to test the 
consistency of the "CMB WMAP concordance cosmol- 
ogy" with the properties of clusters in the low-z Universe. 
With its mosaic of 10 ks overlapping XMM pointings, 
the XMM Large Scale Structure survey (XMM-LSS) has 
been designed to detect a significant fraction of the galaxy 
cluster population out to a redshift of unity over an area 
of several tens of square degrees, so as to constitute a 
sample suitable to address the questions outlined above 
(Pierre et alJl2004|) . The trade-off in the survey design 
was depth versus coverage, keeping within reasonable to- 
tal observing times. 

This configuration allows investigation of the cluster pop- 
ulation down to a depth of about 10~ 14 ergcm~ 2 s _1 
which is comparabl e to the de epest ROSAT serendipi- 
tous surveys (Rosati et al. 2002). However, observations 
are performed with a narrower PSF (FWHM ~ 6" for 
XMM vs - 20 " for the ROSAT PSPC) and very different 
different instrumental characteristics such as background 
noise and focal plane configuration. This led us to de- 
velop a new X-ray pipeline which is presented in the next 
section, along with the principles of the computation of 
the survey selection function. Spectroscopically identi- 
fied clusters undergo detailed measurements, which en- 
abled us to track the evolution of the low-mass end of the 
cluster population (Sec. 3). Source statistics and the L-T 
relation for the CFHTLS Dl area are presented in Sec. 4. 
In the following, we assume a ACDM cosmology and 
give all fluxes in the [0.5-2] keV band. 
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2. XAMIN- A NEW X-RAY PIPELINE 
2.1. Design 

The two major requirements of the XMM-LSS X-ray pro- 
cessing were to reach the sensitivity limit of the data in a 
statistically tractable manner in terms of cluster detection 
efficiency, and hence to calculate the selection function 
of the detected objects. The package that we have de- 
veloped, Xamin, combines the sensitivity of the multi- 
resolution wavelet analysis for source detection with the 
rigour of a likelihood analysis for assessing the signif- 
icance of the detected sources, handling the complex 
XMM instrumental characteristics. Both steps use Pois- 
son statistics. The whole procedure has been validated 
by means of ext ensive simulations of point-like and ex- 
tended sources JPacaud et alJ2005l) . 



2.2. Source classification 

Simulations of the LogN-LogS X-ray point source pop- 
ulation give for Xamin a 90% completeness level down 
to a flux of 4 x 10 -15 ergcm~ 2 s _1 for 10 ks exposures. 
As a rule of thumb, the corresponding sensitivity for 
'typical' cluster sources is 2-3 times higher. However, 
cluster detection efficiency depends not only on the 
object flux and size, and the instrumental PSF, but also 
very much on the background level and on the detector 
topology such as CCD gaps and vignetting, as well on the 
ability of the pipeline to separate close pairs of pointlike 
sources. We thus stress that the concept of sky coverage, 
i.e. the fraction of the survey area covered at a given flux 
limit, is strictly valid only for point-sources because, for 
faint extended objects, the detection efficiency is surface 
brightness limited (rather than flux limited). Moreover, 
since the faint end of the cluster luminosity function 
is poorly characterised at z > 0, it is not possible to 
estimate a posteriori what fraction of groups remain un- 
detected, unless a cosmological model is assumed, along 
with a thorough modelling of the cluster population out 
to high redshift. 

Consequently, with the goal of constructing deep 
controlled samples suitable for cosmology we de- 
fine, rather than flux limits, two classes of extended 
sources corresponding to specific levels of contami- 
nation and completeness. The selection is performed 
in the Xamin output parameter space defined by 
extent likelihood, extent, detection 
likelihood. Selection criteria have been established 
by means of extensive simulations for cluster apparent 
core-radii (extent) ranging from 10" to 100", and total 
number counts from 50 to 1000. The cluster surface 
brightness distribution was assumed to follow a /3-model, 
with j3 — 2/3. The CI class is defined to have "no 
contamination", i.e. no point sources missclassified as 
extended. For the C2 class, selection criteria are relaxed 
such as to allow for 50% contamination by spurious 
extended sources. The classification has been in turned 
checked against some 60 spectroscopically XMM-LSS 
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Figure 1. Cluster selection criteria in the Xamin param- 
eter space populated by simulations. Point sources are 
represented by green diamonds, clusters by blue squares; 
red triangles indicate spurious sources. The horizontal 
and vertical solid lines delineate the C1+C2 class. The 
dotted vertical line indicates the restriction to the CI 
class to which the detection likelihood > 32 
condition was added, guaranteing an uncontaminated 
sub-sample of clusters 



clusters confirmed to date. CI clusters are high surface 
brightness extended objects; this selection inevitably 
retains a few nearby galaxies, but these are readily dis- 
carded from the sample by inspection of optical overlays. 
The C2 sample includes fainter clusters than CI, and also 
a number of nearby galaxies, saturated point sources and 
unresolved pairs, as well as cases badly effected by CCD 
gaps; the contamination is a posteriori removed by the 
visual inspection of the optical overlays as well as by the 
outcome of the spectroscopic identification programme. 
The principle of the procedure is illustrated in Fig. ^ 



2.3. The survey selection function 

In parallel, simulations provide the necessary basis for 
the computation of the selection function. They allow us 
to derive detection probabilities as a function of source 
core-radius and countrate for any exposure time, back- 
ground level and position on the detector (Fig. [2}- As an 
example, we show the dn jdz prediction for the C 1 cluster 
populati on, assuming the fo llowing model: ACDM cos- 
mology (Bennett et al. 20 031) + P(k) power s pectrum with 
transfer function from Bardeen et alj (|19 86[) and comov- 
ing halo number density from me lsheth&^)rmen| (|1999|) 
mass function, il/200 — T relatio n from lArn aud et all 
( I2005I) . L — T relation from lArnaud & Evrard i 19991) 
and a constant core radius of 180 kpc. The cluster/group 
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Figure 2. Detection probability for CI -type clusters as a 
function of count rate and core-radius, averaged for the 
20 inner arcmin of the XMM field, and for an exposure 
time of 10 ks 



population was simulated down to T = 1 keV (no evo- 
lution of the cluster M-T-L scaling laws was assumed, 
as it is currently unconstrained by observations for the 
1 < T < 3 keV range which constitute the bulk of our 
population ). The predicted CI dn/dz is shown in Fig. 
|3J along with the observed redshift distribution of the ob- 
served CI population. The agreement is very satisfac- 
tory and the data suggest a deficit of clusters around a 
redshift of 0.5, probably induced by a cosmic void. Ob- 
served cluster numbers from the currently available 51 
XMM-LSS pointings are 7/deg 2 and 12/deg 2 for the CI 
and C2 class respectively. Finally, we have investigated 
for our simple cosmological model, to what extent the 
CI and C2 classes are comparable to flux limited sam- 
ples (Fig. 0). The main result is that the CI sample ap- 
proaches, in terms of number density, a flux limited sam- 
ple of about 2 x 10~ 14 ergcm _2 s _1 . But (1) more high-z 
objects are detected, while nearby low-surface brightness 
groups are not retained, and (2) the CI sample is strictly 
defined from X-ray criteria (and hence is not contami- 
nated) while it would be necessary to examine a large 
number of sources to clean a putative flux limited sample 



at 2 x 10 1 erg cm 2 s , with many of them not being 
unambiguously characterized as extended or point-like. 



3. CLUSTER MEASUREMENTS 



Each spectroscopically confirmed cluster undergoes de- 
tailed spectral and spatial fits in order to determine its 
temperature, flux and bolometric luminosity within i?5oo 
(standard overdensity radius). In particular, we have 
demonstrated that, under specific statistical model and 
binning conditions, we reach a 20% accuracy in temper- 
ature measurements with ~ 100 a nd ~ 300 photons for 
1 and 2 keV groups respectively (Willi set alJl2005|) . It 
turns out that the CI cluster sample is almost identical to 
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Figure 3. Redshift distribution of the CI clusters for the 
current XMM-LSS area: 48 pointings covering 4.6 effec- 
tive deg 2 ( only the inner 20 arcmin are considered) - that 
is 29 objects. The dotted line corresponds to the predic- 
tions from a simple halo model in a ACDM cosmology 
( see text). 
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Figure 4. Comparison of the CI and C1+C2 redshift dis- 
tribution with various flux-limited samples, assuming the 
simple cosmological model described in the text 



Figure 5. Example of a CI cluster at z = 0.26 (F ~ 
8 x 10~ 14 ergcm _2 s _1 J. The cluster X-ray contours 
are overlaid on a u,r,z CFHTLS composite. They are 
drawn from the co-added [0.5-2] keV M0S1 +M0S2+pn 
mosaic, filtered in the wavelet space using a signifi- 
cance threshold of 10~ 3 for Poisson statistics (not cor- 
rected for vignetting). The intensity scale is logarithmic 
(counts/pixel/second, not corrected for vignetting). Boxes 
indicate cluster members for which we obtained a redshift 
using LDSS at the Magellan telescope. 



the sample for which we can measure temperatures. 



4. RESULTS FROM THE Dl CFHTLS AREA 



The Dl CFHTLS area covers ld eg 2 and constitutes the 
central part of the XMM-LSS (see lPierre etail j2004l) for 
a general layout of multi-wavelength coverage associated 
with the XMM-LSS). It includ es among ot hers the Vi- 
mos VLT Deep Survey (VVDS. lBondi etal] J^OOl l. We 
present a summary of the properties of the cluster sample 
for this region 1 . We have spectroscopically identified 13 
clusters over the 0.8 deg 2 effective area. Out of these, 
8 are CI clusters, 1 is classified as C2, and 4 are clus- 
ters not entering the classification (very faint objects or 
clusters contaminated by a point source). The selection 
function for these 4 latter systems is unknown, but they 
are interesting objects indicative of our ultimate detection 
limit, i.e. ~3-5x 10" 15 ergcm _2 s _1 . Two examples 
are illustrated on Figs. [5] and [6] The C1+C2 clusters 
span the 0.05 < z < 1.05 redshift range, with bolometric 
luminosities ranging from 0.1 to 4 x 10 44 erg s _1 . 
7 clusters are found above a flux of 2 x 10 -14 



'XMM exposure time is 20 ks over the region but we keep the same 
class definition as the S/N increase is only \/2; the net effect being only 
a slight increase in C1/C2 number density ratio for this sub-region 



Figure 6. Example of a faint cluster at z = 0.92 at the de- 
tection limit (F~ 3 x 10~ 15 ergcm~ 2 s -1 ). The source 
was not spectroscopically followed-up. Redshift confir- 
mation comes from the WDS data alone. Same symbols 
and contours as in Fig. \5\ 



erg cm s . This translates to about 8.5/deg 2 . We have 
investigated the L-T relation for the 9 clusters for which 
it was possible to measure a temperature. Results are dis- 
played on Fig. 0as a function of redshift. Malmquist 
bias is here obvious since a temperature was derived only 
for the apparently brightest clusters. The overall corre- 
lation appears to be quite tight and well summarizes the 
ability of rather shallow XMM survey-type observations 
to provide important new insights into the cluster pop- 
ulation. This is particularly relevant for groups out to 
z ~ 0.4 that the XMM-LSS is systematically unveiling 
for the first time. Our intermediate redshift subsample 
of groups (0.9-2 keV for z me d = 0.25) contains objects 
more luminous than predicted by the self similar evolu- 
tion model. According to this hypothesis, the luminosity 
scales as the Hubble constant if it is integrated within a 
radius corresponding to a fixe d ratio with respect to the 
critical density of the universe tVoit 2005). From the lo- 
cal universe, we know that low temperature groups show 
a larg er dispersion in the L-T r elation than massive clus- 
ters (|Helsdon & Ponman 2003). This reflects their indi- 
vidual formation histories, since they are particularly af- 
fected by non-gravitational effects, as well as the possible 
contributions from their member galaxies. The apparent 
biasing toward more luminous objects and/or cooler sys- 
tem could come from the fact that we detect more eas- 
ily objects having a central cusp, i.e. putative cool-core 
groups. This is under investigation using the full sample 
available from the current XMM-LSS area (Pacaud et al 
in preparation). Detaile d presentation of the D 1 catalogue 
and results are given by Pierr e et alJ (120051) . 
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Figure 7. L(R§qo)-T relation for the Dl area clusters 
for which a temperature was obtained. The solid lines 
gives the mean local L-T relation, while the dotted line 
is the expected luminosity enhancement from the self- 
similar model, at a mean redshift of 0.5. The various sym- 
bols show the 3 distinct redshift (and luminosity/mass) 
regimes of the measured clusters. 



5. SUMMARY AND CONCLUSIONS 



We have developed an X-ray pipeline enabling the con- 
struction, from shallow XMM pointings, of samples of 
galaxy clusters with controlled selection effects. In par- 
ticular, "first class" clusters constitute a sample selected 
upon X-ray criteria alone (once nearby galaxies are re- 
moved). This approach offers the advantage of avoid- 
ing hypotheses about the faint end of the cluster LogN- 
LogS (currently not explored beyond the local universe) 
since this can become critical for faint samples said to be 
"flux limited". The CI class contains objects as faint as 
~ 1CP 14 erg cm _2 s _1 out to redshifts of unity. The final 
C1+C2 sample reaches a density of 12/deg 2 but the initial 
C2 selection has to be cleaned (by inspection of overlays 
and possibly optical spectroscopy) of a similar density of 
spurious cluster candidates (^5/deg 2 ). 
We systematically unveil the low end (T < 2 keV) of 
the cluster population out to a redshift of 0.4. For the Dl 
CFHTLS area, we find a cluster density of 8.5 clusters per 
deg 2 having a flux larger than 2 x 10~ 14 ergcm~ 2 s _1 . 
This is higher than t he 4-5 clusters /d eg 2 given by the 
RDCS LogN-LogS dRosati et aljfl9 98) and the shallow 
XMM/2dF survey llGa ga et al]l2005t) . which finds 7/2.3 
= 3 /deg 2 for the same flux limit. Given the size of the 
studied area, our results are certainly subject to cosmic 
variance which may also affect the 2.3 deg 2 XMM/2df 
survey and, to a lesser extent, the RDCS, which cover s 
5 deg 2 at 2 x 10~ 14 ergcm" 2 s _1 jRosati et al]l2002h . 
For comparison, our simple cosmological model predicts 
some 7.5 clusters /deg 2 hav ing T > 1 keV ab ove a flux 
of 2 x 10" 14 ergcm- 2 s _1 dPacaud et all2005l) . 
We present the first L-T relation from XMM survey- 
type observations over a contiguous area. The rela- 
tion contains 9 clusters out to a redshift of unity over 



only ~ ldeg 2 and so, opens remarkable perspectives 
for the study the evolution of the scaling laws for the 
cluster/group population with moderate XMM exposure 
times. On-going work aims at characterizing the evolu- 
tion of the L-T relation for the T < 2 keV population 
over the currently available 7 deg 2 of the XMM-LSS (in- 
cluding the Subaru Deep Survey). In parallel, extensive 
simulations of XMM images enable us to investigate the 
effect of the cluster scaling laws (M-L-T-R), and of their 
evolution, on the survey selection function. These will be 
further constrained by the APEX and AMiBA Sunyaev- 
Zel'dovich surveys of the region to be performed in 2006. 
It will then be possible to model in a self-consistent way 
the cluster population in parallel with constraining cos- 
mological parameters. 
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